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Abstract. In this special issue article, based on the talk with the same title in 
session B5 (Theoretical and Mathematical Cosmology) at GR19, we review the case 
of vorticity generation in cosmology using cosmological perturbation theory. We show 
that, while at linear order the vorticity evolution equation has no source term in 
the absence of anisotropic stress, at second order vorticity is sourced by gradients in 
entropy and energy density perturbations. We then present some estimates for the 
magnitude and scale dependence of the vorticity power spectrum using simple input 
power spectra for the energy density and entropy perturbations. Finally, we close with 
possible directions for future work followed by some hints toward the observational 
importance of the vorticity so generated, and the possibility of primordial magnetic 
field generation. 
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1. Introduction 

Vorticity is a common phenomenon in situations involving fluids in the 'real world' (see 
e.g. Rcfs. [1, 2J). There has also been some interest recently in studying vorticity in 
astrophysical scenarios, including the inter galactic medium [31 H] , but relatively little at- 
tention has been paid to the role that vorticity plays in cosmology and the early universe. 

In classical fluid dynamics the evolution of an inviscid fluid in the absence of body 
forces is governed by the Euler, or momentum, equation 

|E + ( V . v)« = ~ vp, (i.i) 

at p 

where v is the velocity vector, p the energy density and P the pressure of the fluid. The 
vorticity, w, is a vector field and is defined as 

uj = Vxv, (1.2) 

and can be thought of as the circulation per unit area at a point in the fluid flow. An 
evolution equation for the vorticity can be obtained by taking the curl of Eq. ( II. ip . 
which gives 

du 1 

- = Vx(i;xw) + — Vp x VP . (1.3) 
at p z 

The second term on the right hand side of Eq. (11.31) . often called the baroclinic term 

in the literature, then acts as a source for the vorticity. Evidently, this term vanishes 

if lines of constant energy and pressure are parallel, or if the energy density or pressure 

are constant. A special class of fluid for which the former is true is a barotropic fluid, 

defined such that the equation of state is a function of the energy density only, i.e. 

P = P(p), and so (l/p 2 )Vp x VP = 0. 

For a barotropic fluid, the vorticity evolution equation, Eq. (II. 3p . can then be 

written, by using vector calculus identities, as 

Duj d(jj 

~Dt = ~dt + (V ' V)W = ^ ' V)V ~ ' V) ' (L4) 
which makes it clear that, in this case, the vorticity vector has no source, and so uj = 
is a solution to Eq. (II. 4p . 

Allowing for a more general perfect fluid with an equation of state depending not 
only on the energy density, but of the form P = P(S, p) will mean that, in general, the 
baroclinic term is no longer vanishing, and so acts as a source for the evolution of the 
vorticity. This is Crocco's theorem [5] which states that vorticity generation is sourced 
by gradients in entropy in classical fluid dynamics. 
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2. Elements of Cosmological Perturbation Theory 

In order to study vorticity in the early universe, we must use general relativity. 
Since Einstein's equations are notoriously difficult to solve, and relatively few exact, 
inhomogeneous solutions exist, we invoke the technique of cosmological perturbation 
theory. This means starting with a homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) spacetime as a background and adding on top layers of inhomogeneous 
perturbations. Of course, perturbing the geometry of the spacetime invokes 
perturbations in its matter content, so, for example, the energy density is perturbed 
as 

p{x\n) =p (7}) + 6p(x i ,rj), (2.1) 

where r\ is conformal time and a subscript zero denotes the homogeneous background 
part of the quantity. The perturbations are then expanded order by order in a series as 

$P = $Pi + \^P2 + ■■■ , (2.2) 

and p 3> 5 Pi > &P2i where the subscripts denote the order of the perturbation. In 
analogy with the classical case, we want to consider perturbations of a perfect fluid with 
equation of state P = P(S,p). Expanding the pressure perturbation in a Taylor series 
gives 

dP dP 

SP -as ss + ^"' < 2 ' 3 > 

which can also be written as 

5P = 5P nad + c 2 s 5p, (2.4) 

where we have introduced the non-adiabatic pressure perturbation, <5P na d and the 
adiabatic sound speed squared, c 2 . This can then be extended beyond linear order; 
see Ref. [5] for details. 

The most general perturbations to the FRW line element are given by [7J |8] 

ds 2 = a 2 (r]) [-(1 + 2<p)dri 2 + 2B i dx i dr i + ((% + 2C ij )dx i dx j ] , (2.5) 

where the perturbations of the spatial components of the metric can be further 
decomposed into scalars, vectors and tensors (characterised by their transformation 
behaviour on spatial hypersurfaces) as 

Bi = B )t - Si , = -tpSij + E tij + F(ij) + ^hij , (2.6) 

where Si and F^ are divergence-free vectors, and hij is a divergence- free, traceless rank 
two tensor. Each metric perturbation is then expanded order by order in an analogous 
way to the energy density, Eq. (12. 2p . 
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Splitting the spacetime into a background and perturbation introduces spurious 
coordinate artefacts, or gauge modes, which have the potential to cause a problem. In 
order to remove these gauge modes, we need to make a gauge choice. Since there are 
myriad articles in the literature dedicated to this issue, we do not further comment on 
this here, pointing the interested reader to, for example, Refs. [Zl El El [10] for more 
information. For the rest of this article, we will work in the uniform curvature gauge, 
for which E = = ip and Fi = 0, and neglect tensor perturbations, since they are 
negligible compared to scalar perturbations. This gives the line element 

ds 2 = a 2 (r]) [-(1 + 2<p)dn 2 + 2B i dx i dri + S ij dx i dx i ] , (2.7) 

and the fluid four-velocity, 

u M = ~a [2(1 + 00 + 2 - <f>\ + v lkVl k , -2V U - V 2l + 20 1 B li ] , (2.8) 

where where v % is the spatial three- velocity and V 1 = v l + B l . 

Evolution and constraint equations are obtained, in general relativity, through 
the covariant conservation of energy-momentum, V M T M „ = 0, and the Einstein field 
equations, G^ u = 8nGT^ v , where T^ u is the energy- momentum tensor, and G^ u the 
Einstein tensor, as usual. In cosmological perturbation theory, we obtain such equations 
order by order. As an example, consider the momentum conservation equation, from 
V^T^j = 0, at linear order 

= 0, (2.9) 

where % is the Hubble parameter in conformal time, defined as H = a' /a, the prime 
denoting a derivative with respect to conformal time. Eq. ( 12. 9 j) is the analogue of the 
Euler equation ( II .ip in cosmological perturbation theory. 



v( l + n(i-3c 2 )v ll + 



SPi 
Po + Po 



3. Vorticity in Cosmology 

In general relativity, the vorticity tensor is defined as the projected anti-symmetrised 
covariant derivative of the fluid four velocity, that is [S] 

= V«Vj>u {a . A , (3.1) 

where V^ v is the projection tensor into the instantaneous fluid rest space and is given 
by 

= + u p. u v (3-2) 

Note that, in analogy with the classical case, it is possible to define a vorticity vector 
as uj^ = ^e^yU^, where e^ ul is the covariant permutation tensor in the fluid rest space 
(see Ref. PUDS]). 
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The vorticity tensor can then be decomposed in the usual way, up to second order 
in perturbation theory, as Wjj = Uuj + h^nj- Working in the uniform curvature gauge, 
and considering only scalar and vector perturbations, we can obtain the components of 
the vorticity tensor by substituting the expressions for the fluid four velocity, Eq. (12 .8p . 
along with the metric tensor into Eq. ( 13. ip . At first order this gives us 

ojuj = aVi[ij] , (3.3) 

and at second order 



u 2ij = aV 2 [i,j] + 2a V^V^ + 4> lt [i (Vi + B\)^ - </>iPi[ij] 



(3.4) 



3.1. Vorticity Evolution 

In order to obtain an evolution equation for the vorticity at first order, we take the time 
derivative of Eq. ( I3.3p . Simplifying, and using the first order equations, gives 



h3\ 



-H{1 - Sc^ujuj 
a 



(3.5) 



and so, from Eq. ( 13. 3p . 



, Uj - 3Hc s u Uj = . (3.6) 

This reproduces the well known result that, during radiation domination, |u;iiju;f| oc a~ 2 
in the absence of an anisotropic stress term [5|. 



At second order things are somewhat more complicated. We now need to take the 
time derivative of Eq. (13. 4p . and use the first order evolution and constraint equations, 
as well as the second order conservation equations in order to eliminate all the metric 
perturbation variables. We omit the derivation and the governing equations in this 
review article (see Ref. [13] for details of the derivation and the governing equations), 
but calculations give 

'5Pi + 5pi x 



2ij 



mci^a + 2 



Po 



Pn 



+ vf k 



(3.7) 



+ 2 (V* - X*) u lijih - 2 (X* d - V&) u lik + 2 {X k l t - V^) u ljk 

(Spi tj 5P nadl>i - 5p lti 5P nadld ) 



Po 



Po + Po 

where Xu is given entirely in terms of matter perturbations as 

'A-kGo 2 



H 



(Spiis-nipo + pfiVu) 



{3.1 
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In fact, even assuming zero first order vorticity, that is louj = 0, the second order 
vorticity evolves as 

3Uc s u)2ij = — < SHV 1{i 6P n&dlJ] + l —— > , (3.9) 

Po + -fO L ft + r J 

and so we see that there is a non zero source term for the vorticity at second order 
in perturbation theory which is, in analogy with classical fluid dynamics, made up 
of gradients in entropy and density perturbations. Note that, in the absence of a 
non-adiabatic pressure perturbation, we recover the result of Ref. [11] that there is 
no vorticity generation. 

3.2. First Estimates: Magnitude and Scale Dependence 

In this section we present an estimate of the magnitude and scale dependence of the 
power spectrum for the vorticity, assuming simple power-law input power spectra. The 
results presented here are calculated in detail in Ref. [T2] . 

Working now in the radiation era, and neglecting the vector perturbations (the first 
term on the right hand side of Eq. (13. 9p ). we obtain 

9a 

U 2ij ~ H^ij = — 2$Pl ^-fnadlji] = ^ij j (3.10) 

°Po 

and we can define the power spectrum of the vorticity in the usual way, in Fourier space, 
as 

(uZik^ubfatV)) = ^&{ki-k2)V u {k,ri). (3.11) 
We take the following input power spectra 

= A ( * ) (^)" 4 , SP^(k,v) = D (i) 2 g)" 5 ,(3.12) 

where the time evolution of the energy density is obtained by solving the first order 
equations, and the wavenumber dependence from relating an initial ansatz to Wmap 
data. The non-adiabatic pressure input spectrum comes from demanding that it decays 
faster than the energy density, and that its spectrum is bluer than that of the energy 
density. The amplitudes can also be related to Wmap parameters. These input power 
spectra allow us to compute the power spectrum for the vorticity analytically, giving 

V (kn)- —k 5 ^\n 2 a{ko) Vfe- 12 A 8 

I wvAj V) ~ 256 fc o 47r4 In ^ _ a(ko) ) ^ 

16 kt ( k \ 7 12 kl ( k \ 9 4 / k \ 11 ' 

(3.13) 



x 



135^ \k J 245 k 2 Q \k J 1575 V.V 
where k c denotes the large wavenumber or small scale integration cutoff, and the other 
parameters are given in Ref. [TJ]. Substituting in an illustrative value for the cutoff, 
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V UJ (k,r ] )= V 2 \n 2 (^- 



Q-20 / * 



k c = lOMpc" 1 , and using the values of the Wmap parameters from Ref. [H], taking a 
conservative value of 10% of the maximum value for a, we obtain the following vorticity 
power spectrum 

0.87 x 10~ 2 f -M + 3.73 x 10~ n ( ^~ 
k J \k 

n - 

j Mpc 2 . (3.14) 
ko/ J 

This spectrum has a non-negligible magnitude which depends upon the small scale 
cut-off k c and the chosen parameters. As this is a second order effect the magnitude is 
somewhat surprising. As can be seen, the result has a dependence on the wavenumber 
to the power of at least seven for our choice of non-adiabatic pressure input spectrum. 

4. Future Directions 

In the above section we have presented estimates for the vorticity power spectrum based 
on simple, power-law input power spectra. While this is a good first approximation, in 
order to obtain more realistic estimates of the magnitude of the early universe vorticity, 
we need to go beyond the simple ansatz for the non-adiabatic pressure perturbation 
input spectrum. 

One way in which a non-adiabatic pressure perturbation can be generated is through 
the relative entropy perturbation between two or more fluids or scalar fields. For 
example, the relative entropy or isocurvature perturbation, at first order, between two 
fluids denoted with subscripts A and B is 



s AB = m( 5 -^- 5 -^]. (4.i) 

\P0B POA/ 

In a system consisting of multiple fluids, the non-adiabatic pressure perturbation is split 
as 0OTO] 

SP nad = 5P intr + SP ieU (4.2) 

where the first term is the contribution from the intrinsic entropy perturbation of each 
fluid, and the second term is due to the relative entropy perturbation between each 
fluid, Sab, and is defined as 

6Pld ~ 7^T7 Yl PoaPob ( c b ~ c a) s ab , (4.3) 

where c\ and c\ are the adiabatic sound speed of each fluid. Thus, for a multiple 
fluid system, even when the intrinsic entropy perturbation is zero for each fluid, there 
is a non-vanishing overall non-adiabatic pressure perturbation. This can be extended 
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to the case of scalar fields by using standard techniques of treating the fields as fluids 
(see e.g. Ref. [B]). Therefore a possible next step will involve using the relative entropy 
spectrum calculated from multi-field inflation as an input for the non-adiabatic pressure 
perturbation. This will enable us to obtain a more realistic description of induced 
vorticity in the early universe. 

5. Discussion and Conclusions 

In this brief article, we have reviewed current progress in the generation of vorticity in 
the early universe through non-linear cosmological perturbations. We started out by 
considering the familiar case of classical fluid dynamics and showed that the evolution 
equation for the classical vorticity contains a source term which is only non-zero if the 
equation of state is a function of two variables: the energy density and the entropy. We 
then briefly introduced some elements of cosmological perturbation theory and, working 
in the uniform curvature gauge, derived evolution equations for the vorticity. At linear 
order we reproduce the well known result that vorticity decays with the expansion of 
the universe in the absence of anisotropic stress. However we showed that vorticity is 
sourced is at second order by gradients in energy density and non-adiabatic pressure 
perturbations. This is in analogy with the classical case. 

We then presented some first estimates of the power spectrum of the induced vor- 
ticity by using simple power laws as the input power spectra for the energy density 
and the entropy perturbations. The results show that the magnitude of the vorticity 
power spectrum under this approximation is non-negligible and the amplification due to 
the large power of k is huge. Finally, we briefly touched on the non-adiabatic pressure 
perturbation created by relative entropy perturbations in multi-fluid systems which will 
allow us to go beyond the simple approximation in the future. 

A non-zero vorticity at second order in perturbation theory has important 
consequences for the generation of magnetic fields, as it has been long known that 
vorticity and magnetic fields are closely related (see Refs. [16], [T7]). Previous works 
either used momentum exchange between multiple fluids to generate vorticity, as in 
Refs. [HI UHl EH ED E2J [23l El] , or used intermediate steps to first generate vorticity for 
example by using shock fronts as in Ref. [25] . However, we do not require such additional 
steps. Therefore, an important extension to the work presented in this review is to 
consider the magnetic fields generated by our mechanism which could be an important 
step in answering the unknown question regarding the origin of the primordial magnetic 
field. 
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